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lO , 1- Introduction 

T-H ! 

C^^ ' The class of Gibbs point processes is interesting because it allows us to introduce 

^^ , and study interactions between points through the modelling of an associated en- 

■"sj" ■ ergy function. Historical aspects of the mathematical theory are covered briefly 

^^ , in Kallcnberg (1983). When the energy function is parametrized, one among 

many other methods of estimation, is the maximization of the pseudolikelihood. 
Baddeley and Turner (2000) dealt with some practical aspects of such a para- 
metric method and gave a survey on asymptotic results. They noticed that 
some classical examples (such as the area-interaction model, the Multi-Strauss 
}_( ' model, the 2-type Strauss model) do not satisfy the assumptions of the existing 

5^ I asymptotic normality results. This paper aims at filling this gap. 

Many proposals tried to estimate the energy function from the available point 
pattern data generated by some marked Gibbs point processes. If the energy be- 
longs to a parametric family model, the most well-known methodology is the 
use of the likelihood function, see e.g. Mollcr and Waggepctcrscn (2003) and 
the references therein. The main drawback of this approach is that the likeli- 
hood function contains an unknown scaling factor whose value depends on the 
parameters and which is difficult to calculate. An alternative approach relies on 
the use of the pseudolikelihood. This idea originated from Besag (1974) in the 
study of lattice processes. Besag et al. (1982) further considered this method 
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for pairwisc interaction point processes, while Jensen and M0llcr (1991) gener- 
alized it to the general class of marked Gibbs point processes. A general review 
of the problem of statistical inference on spatial point processes including the 
Takacs-Fiksel method (a parametric method based on a characteristic prop- 
erty of marked Gibbs point processes using Palm measure) and non-parametric 
methods can be found in the recent monograph of Moller and Waggepcterscn 
(2003). 

In order to underline our theoretical contributions, let us present the differ- 
ent papers discussing asymptotic properties of the maximum pseudolikclihood 
estimator. The first work was done by .Jensen and Mollcr (1991). The authors 
obtained consistency for exponential family models of marked point processes. 
They mainly considered inhibition and hard-core models. They notably applied 
their results on the marked Strauss process. .Jensen and Kiinsch (1994) and Masc 
(1995) focused on specific models with two parameters -the chemical potential 
and the inverse temperature- which can be viewed as particular exponential 
family models. .Jensen and Kiinsch (1994) obtained an asymptotic normality re- 
sult by first assuming the inhibition or hard-core property and then the finite 
range property. Masc (1995) established consistency for the class of superstable 
and lower regular potentials introduced by Rucllc (1970). Masc (2000) extended 
his work to the context of marked point processes and provided asymptotic nor- 
mality by adding the assumption of finite range. Our goal is to deal with most 
classical models (see Baddeley and Turner (2000), M0ller and Waggepctersen 
(2003) and Bertin et al. (1999b)) that could be interesting for practical pur- 
poses. They have been put into three categories according to their validity with 
respect to the previous existing works: 1) Overlap area point process 2) Multi- 
Strauss marked model, Strauss disc type process. 3) area-interaction, Geyer's 
triplet process, fc-nearest-neighbour multi-Strauss marked model. Let us notice 
first that all the examples belong to the exponential family and satisfy the local 
stability and finite range properties. Due to the parametrization proposed by 
Jensen and Kiinsch (1994) and Masc (2000), they only consider the first cate- 
gory. By considering the exponential family, Jensen and Mollcr (1991) include 
a larger class of models. However, the required inhibition or hard-core type 
assumptions are only satisfied for examples 1 and 2. Examples 3 are only lo- 
cally stable. What remains to be established is consistency for examples 3 and 
asymptotic normality for examples 2 and 3. In this paper, a general framework 
is proposed taking into consideration the previous remarks. Results are obtained 
using the general theory on minimum contrast estimators, e.g. Guyon (1995). 

Section 2 introduces some background on marked Gibbs point processes. Our 
models are defined in Section 3. In the same spirit as Bertin et al. (1999a), pro- 
viding existence results of stationary Gibbs states, assumptions on these models 
are expressed in terms of the local energy function. We also describe examples 
of interest of this work. Section 4 presents the pseudolikclihood method and our 
main results requiring two additional assumptions. The first one is an identifi- 
ability condition ensuring strong consistency. The second one is related to the 
definiteness of the asymptotic covariance matrix of the maximum pseudolikeli- 
hood estimator. These two assumptions allow us to derive a practical asymptotic 
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normality result. They are verified for all the considered examples in Section 5. 
It is the belief of the authors that these assumptions are not restrictive since they 
should be true for every well-parametrized model. Proofs have been postponed 
until Section 6. 



2. Background on marked Gibbs point processes 

For the sake of simplicity, the framework of this paper is restricted to two- 
dimensional marked Gibbs point processes. All the results must remain valid in 
the general d-dimensional (d > 1) case. Define S^ the Borel cr-algebra on R^, B^ 
the set of bounded Borel susbsets of R^ and A^ the Lebesgue measure on R^. 
Denote also by M, M and A™ the mark space and its corresponding cr-algebra 
and probability measure. Let S := R^ x M, S := S^ (g) A^ and /^ := A^ (g) A™ 
denote respectively the state space and its corresponding cr-algebra and measure. 

For shortness, let us denote cc™ ~ (a;, m) for any a; G R^ and any mark 
m & M. and |A| := A^(A) for any A e B'^. In addition, |/| designates the 
number of elements of some countable set /, A*^ is the complementary of some 
set A in R^ and || • || is the ^^-norm. Let us define for all i = {11,12) G Z^, 
d > and p > Aj(d) := {2 e R^, d (ij - \) < Zj < d {ij H- i) , j == 1, 2} and 
B{i,p) := {k e Z^ : \k - i\ < p} with \i\ := max{\ii\,\i2\). 

Let il denote the set of so-called configurations -of marked points- (p := 
{a;™'}ig/ where /is a subset of N and {{xi,mi))^^j is a sequence of elements of S. 
In particular, any element ip ^ Q has the following representation p = X)ig/ ^x"^' 
as an integer- valued measure on $ such that for every F ^ B^, p{F) G N, where 
(S^m is the Dirac measure at some element x™ G S. The subset of i7 with elements 
p satisfying \p\ := </?($) < -l-oo is denoted by fJ/. The space f2 is equipped with 

the cr-algebra T generated by the family of sets -^ (/? G O : p{F) = n\ with n G N 
and F E Bl- For every F G B^ and (/3 G fi represented as (p ~ J2i£i ^x'"^ ' '^^^ 
introduces ipp := Sie/ x'"' gf '^a;'"* which can be viewed as the configuration 
of marked points of ip in F. Furthermore, for every A G S^, pA conveniently 
denotes (^axm- 

A marked point process is a rj-valued random variable, denoted by <f>, with 
probability distribution P on (f7, JF). The intensity measure Np of P is defined 
as a measure on B^ such that for any F G Bf' 



b- 



Np{F) 



fp{F)P{dp):^EmF)). 
Jn 



In the stationary case, Np{F) = vp\^{F) where vp is called the intensity of P. 
A marked Gibbs point process is usually defined using a family of local spec- 
ifications with respect to a weight process (often a stationary marked Poisson 
process with distribution Q and intensity \q = 1). Let A be a bounded region 
in R^. For such a process, given some configuration (^a<^ on A"^, the conditional 
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probability on A is of the form, for any F E J-: 
with the partition function 

JUa 

and i?A ~ {^p ^VL : < Zty{tf) < +00} where 

Let us define the subset of all admissible configurations 
17 := 1^ e O : (^ e HAee^^A} 

and denote by fi/ := 17/ n Q,. Whereas the finite energy function V (cp) (for any 
if &Q,f) measures the cost of any configuration, the local energy V {ipW) (for 
any ip,Tp G flf) represents the energy required to add the points of ip in ip: 

V{ip\ip) =V{^Uip)~V{ip). 

Let us notice that when ■)/) is a singleton {cc™}, we denote by a slight abuse 
V {x"^\ip) instead of V {{x"^}\ip). It is well-known that the collection of proba- 
bility kernels (n.A)AGB^ satisfies the set of compatibility and measurability con- 
ditions which define a local specification in the Preston's sense (Preston (1976)). 
The main condition is the consistency: 

HaH^' = Ha for A' c A. 

Notice that some conditions are needed to ensure the existence of a probability 
measure P related to any local energy V and any weight process that satisfies 
the so-called Dobrushin-Lanford-RucUe (D.L.R.) equations: 

P{F\TAc){ip) = nA((/3, F) for P a.e. ip e il for any A G Sg and F G T. 

For the general theory of Gibbs point processes, the reader may refer to Kallenberg 
(1983); Stoyan ct al. (1995) and the references therein. 

For some finite configuration ip (resp. some set G) and for all x ^ R^, ipx 
(resp. Gx) denotes the configuration ip (resp. the set G) translated of x. Finally, 
in this work a non-marked point process can be viewed as a particular case of 
marked point processes with M — {0}. 
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3. Definitions and examples of marked Gibbs models 

The framework of this paper is restricted to stationary marked Gibbs point pro- 
cesses based on an energy function invariant by translation, V (</?; 9), parametrized 
by some £ 0, wliere is some compact set of W . The model is also assumed 
to belong to an exponential family, i.e. 

v{v-e) = e'v{^), (3.1) 

where v{ip) = (ui((/3), . . . ^Vp{Lp)) is the vector of sufficient statistics. The local 
energy is then expressed as 

F(a;'"|^;0) = 6^i;(a;"», (3.2) 

where v {x™'\ip) = (wi(a;™|(p), . . . ,Vp{x™'\ip)) := v{ip U {x"^}) — v{ip). 
Our models satisfy the general condition [Mod] described by the following state- 
ments: 

[Mod:S] Stability of the local energy: there exists K > such that for all 

(m, ip) eMx Qf 

^(O^lvs;©) > ^K. 

[Mod:L] Locality of the local energy: there exists D > such that for all 
{m, if) eMxQf 

y(O'"|^;0) = y(O™|<^B(o,i^);0), 

where B{x,r) denotes the ball centered at a; G M^ with radius r > 0. 
[Mod:I] Integrability condition: for i ~ 1, . . . ,p, there exist kJ*'"'^ > 0, ki € N 
such that for all {m,ip) eWL x Qj 

V^{0"'\v)<^t''^\mO.D)\'''■ 

Let us notice that, unlike [Mod:I], the assumptions [Mod:S] and [Mod:L] 
cannot be directly expressed only in terms of the sufficient statistics. Neverthe- 
less, [Mod] is satisfied as soon as for i = 1, . . . ,p, there exist k^'" , k\^'^^ > 0, 
ki G N such that one of both following assumptions is satisfied for all {m, (p) G 
M X O/: 

[Mod-1] 

0., > and - /.ff) < z;.(0™|^) = vdO"^\ipBi,,D)) < «:f"'V8(o,D)l''. 
[Mod-2] 

Indeed, let /i and I2 be the partition of {1, ... ,p} such that for any i e /i (resp. 
i G I2), Vi satisfies [Mod-1] (rcsp. [Mod-2]) then 

F(O™|^;0) = ^0,z;,(O™|^) + ^a,i;,(O™|^) 
ieii ieh 



> -p Uax(0,«;f"f)) -max (|^,;| x max(Kf"'\ /.f "p')) j := -K{e) 
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which ensures [Mod:S] with K := s\\y>s^@K[0) ([Mod:L] and [Mod:I] are 

clearly satisfied). 

Let us also point out that 

• the well-known characteristics [Mod:S] and [Mod:L] associated with fi- 
nite energies that are translation invariant ensure the existence of station- 
ary measures (see Bertin et al. (1999a)). 

• the local stability implies the Ruelle-bound correlation function (see Ruolle 

(1970)) leading to: E {\'^b(o,d)\ ) < +oo for any /3 > 0. 

• A key-ingredient of our proofs is the following one: for any a > 0, any 
9 £ @ and any i — 1 , . . . , p 



E (|t;,(0^^|$)|"e-^"^(°"l*)) < +00. (3.3) 



This condition is fulfilled under the assumptions [Mod:S] and [Mod:I]. 

Let us now present some examples. Except the one based on the fc-nearest- 
neighbour graph, all examples are really classical and can be found e.g. in 
Baddclcy and Turner (2000) and MoUcr and Waggepetcrsen (2003). For each 
example, we present the model through the sufficient statistics, the set of the 
parameters values (including 0) for which the model is defined in the littcra- 
turc and then we verify that [Mod] is satisfied. This proves in particular the 
existence of stationary Gibbs states in M^. 

First of all, note that when Vi{0™'\(p) := 1, then Vi obviously satisfies [Mod- 
2]. Recall that for a non-marked point process M = {0}. 

Overlap area point process 

• This non-marked process is defined for p = 2 and i? > by 

z;i(^):=|^|andt;2(^):= ^ \B {x,R/2) n B {y,R/2)\ , 

{x,y}eV2iv) 

where Vkif) [k > 1) is the set of all subsets of (p with k elements. Alter- 
natively, 

vi{0\ip) := 1 aiidv2{0\ip) := ^ |S(0, i?/2) n S (x, i?/2)| . 

xdip 

Let us notice that 

\B{0,R/2)nB{x,R/2)\ 

= i (^R^Arcos^^ - ||x||v/i?^HNF) l[o,fl,](l|x||). 

• eK xM+. 

• V2 satisfies [Mod-1] with k^'"^^ = 0, k^™p^ = ^, D = R and k2 = 1. 
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Multi-Strauss marked point process 

• Let M = {1, . . . , M}, A™ the uniform probability measure on M andp > 2. 
Decompose ip = U^^-^ip"^ with ip"^ := {cc"* S fi : x™ S t/?} for any m e M. 
The finite energy of this process is defined by 

M 



Vip;e) := J2 0™^'"^w™^'"^((/j"O 



I<mi<m2<-A/ 4=2 

with for all mi, m2 G M and i ~ 2, . . . ,p 

XI iiD-iv^'.-Dr'''"'!^!^!"^^!!) (3.4) 

where < _C)™i'™2 ^ £)mi,m2 <; <; D"iil^^2 < +00. In particular, the 
vector 9 could be ordered as follows: 

9 = i9\9^, ..., 9^') where 6>™ = (6/"'™, 6/"'™+\ . . . , 6/"'^^) 

with 

if TTLi = 77l2 



imi,Tn2 


' fgmi,m2 
J 


Qnii ,m2 


nmi,m2 




(nmi,m2 


nmi,m2 
. ^3 


/)mi,m2 



9' . ,_ _ _ ^ 

otherwise. 

where p = M + >^ (j/ni,m2 _ -^y q^^^ then derives the expression 

l<mi<ni2<M 

of the local energy 

V (0™1 1^; 9) = 0^''™' + H H 0™l.™2^mi,m2(o™i |^™2) 

m2=l i=2 

where for convenience 0™i-"'2 g^^^j Qrn2,rni ^j^jj^q^-j, ^-j-^^ same parameter. 
For all mi,m2 G M and i = 2, . . . ,p™i'™=: 6'"i''"i e R and 61™^'"" G 
wheni?™^''"^=^>0 
M+ wheni^T^^^"^ ^0. 
• When i^™i>™= = and C''"' > 0> ^'"''"" satisfies [Mod-1] with 

^(inf) _ n (sup) _ -. pj _ 7-)mi,m2 J jL i 

i(mi,m2) ~ "' 4(mi,m2) - J^ , -t^ — max„jj<m2 IVp^^.^j dno ^^(mj^m^) - i, 

where ^™i'™^ = 6i(^jni,m2) is the *('^i, ?n2)^th element of the vector 0. Un- 
der the hard-core assumption £)™i''"2 = j > g, y™'^'"^^ satisfies [Mod-2] 
with K^)" \ = 0, K^?"''-^ , := r-rrl and D = max„,,<™, D™„i:™?2- 
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k— nearest-neighbour multi-Strauss marked point process 

• This marked point process is defined similarly as the multi- Strauss marked 
point process except that the complete graph P2{f) in (3-4) is replaced 
by the fc-nearest-neighbour graph (fc > 1). 

• For all mi, ma G M and i = 2, . . . ,p™i.™2. grm.m^ ^ ^^ 

• In Bcrtin et al. (1999b), it is proved that t;™!'™^ satisfies [Mod-2] with 
K ." N = K.™^' , = 13fc and D = 2 max™, <„,, D"L\' ™ . 

Strauss type disc process 

• Let M = [0, M,„ax] with < A/max < +00, A™ the uniform probability 
measure on M and p = 2. This model is defined by 

vi{if) = \ip\ a.nd V2{ip) = ^ l[o,rm+m2](.\\x2 - Xi\\). 

Alternatively, 

z;2(0"|^)= Y. 1[0.™+™'](II-^II)- 

• eM x]R+. 



• 



V2 satisfies [Mod-1] with Kj ~ 0, K2 = 1, D = il/max and fc2 = 1. 



Geyer's triplet interaction point process 

• This non-marked point process is defined for p ~ 3 and i? > by 

V {ip; 6) = 6li|(p| + eiafaM + e^iV^iip) 

where 

V2{^)= Yl l[o,i?](lki -a:^2||) 

{xi.X2}&'P2(<p) 

and 

^(iViiv) {xi,X2}(^V2(0 

Note that V'i{ip) represents the number of triangles of ip with edges of 
lengths lower than R. Alternatively 



V2 



(0|^) = ^l[o,fl](||x||) 



and 

VMV)^ E n l[aM]{\W-X2\\). 

{x,v}&V2[v) {xi,X2}eV2({x,yfi}) 
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• 6» e ]R2 X E+ \ {0}. 

• When 62 > 0, V2 and V3 satisfy [Mod-1] with D ^ R, k'^"^^ = k^'"^^ = 0, 
^(sup) _ ^^^'"P) — i^ /j2 ~ 1 and ^3 = 2. When 02 < 0, V2 satisfies 
neither [Mod-1] nor [Mod-2]. However, Geycr (1999) proved that the 
local energy is stable and local (i.e. [Mod:S] and [Mod:L]) and [Mod:I] 
is satisfied with D ~ R, Kj "^ = K3 ^ = 1, /c2 = 1 and k^ ~ 2. 

Area interaction point process 

• This model is the one-type marginal of the two-type Widom-Rowlinson 
model. Let p = 2 and R > 

V{^;e) = 9i\^\+92V2{^), with V2{^):=\U.,ev'3{x,R)\. 

Note that V2{f) represents the area of the union of discs of radius R 
centered at the points. Alternatively, 

V2i0\ip) := U^e^gj,_,^ju{o}S(x,i?)\U^e^g(,,^jS(x,i?) 



• 



V2 satisfies [Mod-2] with k^2 = 0, k^2 = ""^^ and D ^ 2R 



Remark 1. Jensen and M0ller (1991) have already proved the consistency prop- 
erty in the inhibition case, i.e. [Mod:S] with K =: Q. In particular, they did not 
consider the area-interaction model with negative parameters 61 and 62 (for in- 
stance). This gap has now been filled by extending the result in the case when 
[Mod:S] is satisfied. However, unlike these authors, we require the additional as- 
sumption [Mod:Ij. In order to simplify our assumptions, we deliberately decided 
not to propose this particular case since this last assumption is not restrictive 
and is satisfied for all examples considered above. 

Remark 2. Note that unlike the Multi-Strauss marked point process, neither in- 
hibition nor hard-core assumption is required for the k— nearest-neighbour multi- 
Strauss marked point process since its local energy is naturally stable. 

Remark 3. Concerning the Geyer's triplet process, the case 63 — is not 
considered since it is a particular case of a multi- Strauss marked point process. 

Remark 4. Through these different examples, one may note that some param- 
eters are assumed to be known: for example, the parameters £)™i'™2 j-^^ ^^g 
multi-Strauss marked point process, the hard-core parameter 6, the parameter 
Mmax for the Strauss type disc process, the parameter R for the Geyer's triplet 
process or the area interaction point process, .... Their estimations could be 
investigated by using ad hoc methods. 
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4. MPLE: presentation and asymptotic results 

4.1- Pseudolikelihood 

As specified in the introduction, the idea of maximum pseudolikelihood is due 
to Besag (1974) who first introduced the concept for Markov random fields in 
order to avoid the normalizing constant. This work was then widely extended 
and Jensen and M0llcr (1991) (Theorem 2.2) obtained a general expression for 
marked Gibbs point processes. With our notation and up to a scalar factor, the 
pseudolikelihood defined for a configuration ip and a domain of observation A is 
denoted by PL\ {(p',0) and given by 

PLA{^;e) =cxp( - f e-^^'^l^'^Vlrf^")) n e^^("'"l^\""'^). (4.1) 

V JaxM J ^™g^^ 

It is more convenient to define (and work with) the log-pscudolikelihood, denoted 
by LPLa (</.;©). 

LPLa ((^; 6>) = - /" e-^(^'"l^'^V(c^a;")- J2 V {x"'\ip\x''';e) . (4.2) 

Our data consist in the realization of a point process with energy func- 
tion V{-;6*) satisfying [Mod]. Thus, 9* is the true parameter to be esti- 
mated and it is assumed that 6* G &. The Gibbs measure will be denoted 
by Pg*. Moreover the point process is assumed to be observed in a domain 
A„ © Z?^ = U2;gA„'B(a;, Z)^) for some D^ > D. For the asymptotic normality 
result, it is also assumed that A„ C M^ can de decomposed into Uig/^Ai where 
/„ = IB (0, n) and for i = {ii, 12) G 1? , A^ = Ai(Z3) for some D > Q fixed from 
now on. As a consequence, as n ^ +00, A„ -^ M? such that |A„| -^ +cxd and 

|A„| ^ ■ 
Define for any configuration (/s, f7„ {v\9) = —■r^LPL\^ (lp;9). The maxi- 
mum pseudolikelihood estimate (MPLE) denoted by 0„((p) is then defined by 

9n{f) ~ argmax LPLa^ (ip; 9) = argmin Un i^', 9) . 
9e& " 9e& 

We will also need the following basic notations: 

• Gradient vector of [/„: Un {^fi', 9) := — |A„|~-'^LPi^ {(p; 9) where for any 
bounded Borel set A, I LPL\^ (ip; 9) 1 is defined for j = 1, . . . , p by 

(lPL^^\p,;9)) ^ f z;,(.T"|vp)e-^(-"l'^^%(dx™)- ^ «,(x"|(^\x™) 
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Hessian matrix of C/„: U)^'{(p: 6) := — \An\^^ LP L j^' {(p;9) where for any 

(2; 

■A 



bounded Borel set A, ( LPL J^ ' (ip; 6) ) is defined for j, fc = 1, ... ,p by 



]M 



^ '' i.k ./ A vM 



J'*: JAxM 

Finally, note that from the decomposition of the observation domain A„, one 
has 

and 

U^\ip; 6) = \A„\-' E ^^A, (^; 0) . 
iein 

4.2. Asymptotic results of the MPLE 

This section provides consistency and asymptotic normality of the maximum 
pseudolikelihood estimator. Let us first consider the following assumption 

[Ident] Identifiability condition: there exists Ai, . . . ,Ai, i > p events of fl and 
Af, . . . , A^ events of M such that: 

• the £ events Bj := A™ x Aj are disjoint and satisfy A™ (g)Pe* (B^) > 

• for all ((mi, ipi), . . . , {mi, ifg)) E Bi x ■ ■ ■ x B^ the {i,p) matrix with 
entries Vj{0™'\(pi) is injective. 

Theorem 1. Under the assumptions [Mod] and [Ident], for Pg* —almost every 
if, the maximum pseudolikelihood estimate Oni'fi) converges towards 6* as n 
tends to infinity. 

For the next result consider 

[SDP] For some A := U^ Bfo [-^1 "l ^*(-^) ^^^^ P > 0, there exists Aq,. . . ,Ai, 

£ > p disjoint events of 51 := <ip £ fl : f^.ijj) = 0, 1 < |i| < 2 £= > such 
that 

. forj = 0,...,£, Pe*(Aj)>0. 

• for all {ipo, . . . , ipi) G Aq X ■ ■ ■ X Ai the {£,p) matrix with entries 
LPL^^^ {^,-e)) - [lPL^^^ (^o;^)) . is injective. 

Theorem 2. Under the assumptions [Mod] and [Ident], we have, for any fixed 
D , the following convergence in distribution as n —^ +00 

|A„|i/2[/(2)($.0*) f0„(a>)-r)^AA(O,S(r)), (4.3) 
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where 

s(r)= Y. ^(^^<!(i)(*;ni^iSi)(*;nn- (4-4) 

iGB(0,r-Dl) ^ ^ 

In addition under the assumption [SDPj 

|A„|i/2s„(<i>;i5^5,0„(cI>))-i/2c/CT($;0„(<i,)) (©„($) -0*]^ A/- (0,/p), 

(4.5) 
where for some 9 and any configuration if, the matrix 'Tj_^{Lp; D^ , D, 9) is defined 
by 

S„(y,; D\D, 6) = |A„ri Y. E ^^^a! (v'; 0) LPlI] (^; ^f . 

(4.6) 

Remark 5. Lei us underline that the scaling that yields to asymptotic normality 
corresponds to the usual parametric rate. Indeed, in the d-dimensional case one 
would obtain 

|A„|i/2 = |A„(5)|i/2 ^ |/„|i/2 X I A,(5)|i/2 ==5''/2((2n+l)^)i/2 ^ {2bY/\n''Y'\ 

Remark 6. We would like to underline that S(0*) = S(i5,0*) = S(1,0*) 
where for all D > 

UD,9*)=D-'~ Y e{lPL^II{<^-9^)LPlI]{<^-9^)\ (4.7) 
ieBi 0, e 



5. Back to examples 

This section is devoted to proving that aU of our examples satisfy both assump- 
tions [Ident] and [SDP]. For the assumption [Ident] V_ denotes the matrix 
with entries Wj(0™' \^i) where (ttt,,;, Lpi) £ Bi have to be defined according to the 
different examples. The assumption [SDP] may be rewritten for all k = I, . . . ,£ 
and for all ipk € ^fc and ipo G Aq: 

(yy e K^/ {lpl^t ivk;e*)^LPL^^'> (^o;r))^'^(i(^fc;r)-ii(<^0) = o)^y = 0, 
where for any configuration ip Cz Q and ipo G ^o 

L{ip;9*) := / v {x"'\ip) e-^*'^''^'''"\^'>ii{dx"') 

J AxM 

V [x'^lifo) e~^*''''(^"l'^«V(rfa;'") 



'AxM 
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Concerning this assumption, wc choose D > D in all our examples. 

5.1. Overlap area point process 

Assumption [IdentJ 
Consider 

Ai := {if en-. ipi3{o,D) = 0} 

A2 :- {ipen:ipB(„,D)^{z}.zeB{iO,D/2),D/4)}. 

We have for all {ipi,ip2) e Ai x A2 

1 
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V = 



1 l'2(0|<y92) 



For every ip2 £ A2 such that iip2)t3{o,D) = {z} with for aU z e S((0, D/2),D/A), 
one remarks that \B {0, R/2) nB {z,'r/2)\ := g2i\\z\\) > 0, then det(F) ^ 0. 

Assumption [SDP] 

Denote by ^0 any configuration set. Consider An{'q) for n > 1 and for some 
< rj < D the following configuration set 

^niv) = \f eft : '^Ao("d) = {^i'-- -7^"} "^ith 21, •■ -i^n e B{0,r])\ . 
For any (^„ e A„(77), we have 

ifi=l 



Vj{x\^n)e-''^ ''(^I'^-^d.T 



r|A|e^_ 
- \n7ri?2|A|e^ if j==2. 



where K conies from the local stability property. Let us also remark that 
y^ R{ipn) = nyi + y2 X ^ V2{x\Lpn\x) - 1^ ^ v{;x\ipo\x) 

with 

< n(n - f)g2(»y) < ^ W2(a;|(^„ \a;) < n(n- 1)- 

Therefore by combining these arguments, for every e > 0, we have for n large 
enough 



7ri?2 



\y2\g2{n) < 



3^72/2 Yl V2{x\lfn\x) 



n{n — 1) 



XeifnHA 



—, -T y2 V V2{x\ipn\x)+lJ^{L{Lpn;e*)~R{ip^))\ 

n[n — 1) I ^-^ _ / 

\ xeip„nA / 



<e. 
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By choosing e 



\y2\g2{ri) 



, this leads to y2 = 0. Then, for every e' > we may 



obtain for n large enough 



m\ 



yi + ^iyi,0f{Liip„;e*)-R{ipr,)) 



<e'. 



By choosing e' ~ \yi\/2, this leads to yi — 0. 



5.2. Multi-Strauss marked type models 

Assumption [Identj 

Define for any m,mi^m2 G {!,..., A/} with m2 > m-i and for any i = 2, 

„mi,m2 

• • • 7 y 



4™ 

Ami,m2 



= {fen-. (fiB{0,D) = 0} 
= {m} 

= {^en-. VBio,D) = {2™n, with z e ^(0,^^''"') \^(o, ^rr^)} 



The following events B"''^'™'^ are defined for any i = 1, . . . ^p™i^™i when nii 
1712 and any i = 2, . . . ^p™i'™2 when mi < m2 such that: 



B 



mi.m-2 _ J ^mi ^ ^0 



if 7711 = "T-2 and i = 1 



A™^ X A^»i'™2 otherwise 



One may order these t = p events as Bi, . . . ,Bi where B^j^^.ra^ ;— B^ 



with 



7712 — 1 7771 — 1 KI 

fcI"-'"^ = (mi-l+<5„„„,)+|i-l+ £ (p"-™^ -1)+ Y. Y. (P"'-™^ -1)| 



The corresponding matrix is then V = 



1^77 



/ 1 ^77,777. 

: 

Vi ... 



\ 



K7,M / 



/ Vi 

'•. 
V ... 



and Vrn 



... \ 



Vm I 

/o .. 

1 ■■ 

■■ 

V .. 



with 



o\ 



■•• 

1 / 
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Assumption [SDP] 

Let us first introduce tlie following sets for any 77 > and d > 



A(77, d) = I (zi, Z2) e Ao(i^)' : zi e Z? (^(0, 0), ^ j and Z2 G S ( (d, 0), ^ 

A-(?7,d) = {(^i,^2)eAo(:D)2:ziGZ?((0,0),|) and Z2 e S ((d - |, 0), | 

C A(r/,d) 
^+('7,0?) = |(zi,Z2)eAo(:D)2:ziGS((0,0),|) and Z2 e S ((d + |, 0), | 

C A(T],d). 

For any i G {2, . . . jp™!'™^}, when 77 is small enough, the couple of points 
(zi,Z2) e A(7?,i?r-"^) (resp. A_(7?,i5r-'"^) and ^+(77,7?™-"^)) are such 
that D,"1V'"^ < Z?™"'"' -^ < '^(^1.^2) < A™"'"' +^ < A"' (resp. 

i^rr^ < A"""' -'^ < ^(^1:^2) < A"""' and i?r"'"' < ^(^1,^2) < 

We now derive the following events 
Ao := jyj e SI : </?Ao("D) = 0} 

^r""^('?) = {^ef^:^^„(U) = K'S^2"^nwith(zi,z2)eA(,7,Dr''"=)} 
<'-"=W = {^en:^^^(^) = {zr,4"nwith(zi,.2)eA_(,7,i5r'"=)} 

<r"('?) = {¥'en:^^^(^) = {zr,^2"nwith(zi,z2)eA+(,7,A'""™^)} 

First of all note that for any (^o G ^Oi 

M 

y'LPLf (^0; r) = Y. yr'^e-^^^-^IAI - 

rn— 1 
For any (^"l'"^ G ^™^,"'2 ^^^ ^^^ ^rni,m2 ^ ^rni,m2 for i ^ 2, . . . ,p™i.'"2 

If Rift^l ) = 2/1 +2/1 +22/i '■ "(l-(S..,p™i,™2) 

We leave the reader to check that for every e > there exists rj > small enough 
such that 
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Therefore for every e > we have for rj smah enough 
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<£. 



By ehoosing e = \y^i^2 \i this leads to 2/!^i™^2 ~ 0- By iterating this argument, 



obtain that for any mi, 7712 G {!,..., Af }, y™ 



yp 



0. It 



remains to prove that y-^^' ==••• = 
configuration set indexed by n > 1 



M.M 



y^ ' =0. For this, consider the foUowing 



AT={^en: <^^^(^) == {z^\ . . . , z™^} with zi, . . . , z„ G Ao(i?)} 



For any (yS™^ G A™^ , we have 

M 






7711 ■'ITT'l 

nyi 



Hence for every e > we have for n large enough by using the local stability 
property 



n 

mi — l 

By choosing £ = l^/^""^'"^^ |/2, this leads to y"^^'^"^ = 0. 



?n 1 , rn 1 

Vi 



5.3. k— nearest-neighbour multi-Strauss marked point process 

Assumption [Ident] (resp. [SDP]) is proven without any change in the proof 
of the multi-Strauss marked point process for every fc > 1 (resp. k > 2). The 
proof of [SDP] for fc = 1 is omitted. 

5.4. Strauss disc process 

Assumption [Ident] 

Consider 

Ai := [ip ^n-. (pB(o.D) =0} 

A2 := {(^Gr!:^e(0,i?)={^™'}, 2G6(0,i?/2), m'e[Z?/2,i?]} 
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and define Bi ■.— 'Mx Ai and ^2 := M x A2. Then, for any (mi, (pi, m2, ip2) G 
Bi XB2 

which is injective since det(Vl) = 1. 

Assumption [SDP] 
Consider for n > 1 

Ao := |(^ e SI : </?Ao("d) = ^] 

zi, . . . , z„ e B{Q, D/2), mi, ... , m,, £ [D/2, Z?] j. 
Note that for every ipo € Aq and any ipn € A„ 

■t/Riipn) = nyi + n(n - 1)2/2 

Note also that from the focal stabihty property \i/^L{ipn',d*) < 2(|yi| + |y2|)|A|e-^. 
Then, for every s: > we have for n large enough 



12/21 = 



— -y^ {Liipn;d*) - RniVn)) + 2/2 

n(77, — 1) 



< £. 



By choosing e = I2/2I/2, this leads to y2 = 0. Then, for every e' > we have for 
n large enough 



\yi\ 



-(2/1, of {L{ip.^; e*) - R„i^„)) + 2/1 
n 



<e'. 



By choosing e' = |2/i|/2, this leads to yi — 0. 

5.5. Geyer's triplet point process 

Assumption [Identj 

By considering 

Ai -.^ {ip e il : (PB(0,D) == 0} , 

A2 -.^ {(f e ft : (PB{0,D) = {z}, Z e S((0, D/2),D/4:), ipi3{0,2D)\BiO.D) 

and 



^3 := {(^ e f^ : (^i3(o,i3) = {^1, Z2}, ^1, ^2 e 6((0, D/2),D/4), fB{o,2D)\B{o,D) = 0}, 
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we have for any ((^i, Lp2, ^3) ^ Ai x A2 'X- A^ 
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V = 



1 

111 

1 2 3 



which is clearly injective since det V^ — 1. 



Assumption [SDP] 

Denote by Aq any configuration set. Consider An{ri) for n > 1 and for some 
< 7] < D the following configuration set 

Aniv) ^ \f eft : fAoiD) = {^i:-- -'^"l "^ith zi,.. .,z„ e S(0,77)| . 
We leave the reader to chek that for j = 1, . . . ,p and for any ipn G An{r]) 



Ja 



<nJ-i|A|e^, 



where K comes from the local stability property. Let us also remark that 



l/R^ipn) = nyi + n{n - 1)2/2 + 



n{n- l)(n-2) 



x^ipoHA 



Therefore by combining these two arguments, for every e > 0, we have for n 
large enough 



m\ 
2 



7/3 lf{Li^n;e*)~Riipn)) 

2 n{n-l){n-2) 



<e. 



By choosing e = ^^, this leads to 2/3 = 0. Then, for every e' , we may obtain for 
n large enough 



12/2 1 



y2 



{yuy2,of{L{^n;0*) - R{^n)) 



nin — 1) 



<e'. 



By choosing e' ~ |j/2|/2, this leads to 1/2 = 0. And then for every e" for n large 
enough 



|yil = 



yi 



{yi,O,oy{L{^n;0*)-Rivn)) 



<e" 



By choosing e" — ^^, this finally leads to yi = 0. 
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5.6. Area-interaction model 

Assumption [IdentJ 

By considering 

Ai -.^ {ip en-. (pB{o,D) = 0} 

and for some small rj > 

A2 :=yl2(r;) = {ifi E Q : ipB(Q.D) = {z}, z e B{0,ij)} . 
we obtain for any (ipi, 1P2) G Ai x A2 

1 7ri?2 



' 1 V2i0\(p2) 

Since for any r; < i?, < V2iO\ip2) < nR'^, det(F) ^ 0. 

Assumption [SDP] 

Consider Ak{rj) for fc = 0, 1, 2 and for some < rj < D the following configura- 
tion set 

Akin) = [v ^^- ^Aa(D) = {zi, . . . , Zk+i} with zi,...,Zk+i e B{0,r])j 

For any ipk e A-(??) for fc = 0, 1, 2, 

i?(V3i) = 2yi + y2.gi(vi)) - (yi + V2T^R^) 
R{ip2) = 3yi + y292{^2) - (vi + V2T^R^) 

For every e > 0, there exists 77 > such that for any i^k G ^fe('?) (^ = 1:2), 
\'ifL{ipk]d*)\ < e and such that \y2gk{Vk)\ < £■ Then, 

\yi\ = \yi - y^iLiipi;e*) - L2((p2; ^*)) - (i?(</'i) - Ri'P2))\ < 4e. 

By choosing e ~ \yi\/8, this leads to yi ~ 0. Now, 

\y2\TTR^ - |y27ri?2 - (0,y2f (^(¥'i;0*) - R{m))\ < 2e. 
And by choosing e = 7ri?^|j/2|/4, this finally leads to 2/2 = 0. 

6. Annex: Proofs 

6.1. Tools 

Let us start by presenting a particular case of the Campbell Theorem combined 
with the Glotzl Theorem that is widely used in our future proofs. 
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Corollary 3. Assume that the (marked) point process $ with probability mea- 
sure P is stationary. Let K be a hounded Borel set, let tp E fl and let g be a 
function satisfying g{x™,ipx) = 5(0™, <p) for all a;™ £ S. Define M a random 
variable with its distribution A™ and f{m,ip) = g{Q"^,i.p)e~^^^ ^'^> and assume 
that f eL^{X^®P). Then, 

Y^ g{x^-, $ \ X™)) = I A| ^( g (0^^ $) g-^Kl*) ) (6.1) 



E 



Proof. 



ai^e^A ^^ K^Gip 



I I g{x^,^)U{x)PU{d^)Np{dx^) 

[ [ g(x™,^)lA(x)j.p(x™)P2„(d^)/i(dx'") 
's Jn 

I f g(.T™,(p)lA(x)e-^(^'"l'^)p(d^)/.(da;'") 
' 5(x"\(^)e-^(=^'"l'^)p(d^)M(da;'") 



AxM Jn 

= |A| / /.g(a;'",(^)e-^(^'"l'^)p(d^)A"(dm) 
Jm Jo 

= |A|^(5(0'^^<l>)e-^Kl*)' 

where !^p(-) is the Radon-Nikodym derivative of Np with respect to /i. ■ 

Let us now present a version of an ergodic theorem obtained by 
Nguyen and Zessin (1979) and widely used in this paper. Let D > and denote 
by Ao the following fixed domain 

Theorem 4 (Nguyen and Zessin (1979)). Let {Hg,G G B^} be a family of 
random variables, which is covariant, that for all x gM.^, 

HgA'Px) = Hg{(p), fora.e.if 

and additive, that is for every disjoint Gi, G2 E B^, 

Hgi\jG2 = -f^Gi + Hg2, a.s. 

Letl he the sub— a— algebra of T consisting of translation invariant (with proba- 
bility 1) sets. Assume there exists a nonnegative and integrahle random variable 
Y such that \Hg\ < Y a.s. for every convex G C Aq. Then, 

^^5* T7^^G„ = T-r-T^(i^Ao|2^), a.s. 
ri^+00 G„ Ao 



for each regular sequence Gn 
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6.2. Proof of Theorem 1 

Due to the decomposition of stationary measures as a mixture of crgodic mea- 
sures (see Preston (1976)), one only needs to prove Theorem 1 by assuming 
that Pe* is ergodic. From now on, Pg* is assumed to be ergodic. The tool used 
to obtain the almost sure convergence is a convergence theorem for minimum 
contrast estimators established by Guyon (1995). 
We proceed in three stages. 
Step 1. Convergence of Un{^',d)- 

Decompose Uni(p; 0) = ^ (-ffi,A„ if) + ^2,a„ (<^)) with 

ffi,A„(^)= / e-^(^'"l'^^^V(da;'") 

J AnXM 

and 

Under the assumption [Mod], one can apply Theorem 4 (Nguyen and Zcssin 
(1979)) to the process Hi,a„- And from Corollary 3, we obtain Pe*— almost 
surely as n ^ +cxd 



ifi,A„ -^eU e-^(°"l*^^'A"(dm)) = £;(e-^Kl*^^)). (6.2) 




Now, let G C Aq, we clearly have 

\H2M^)\< E \V{x"^\^\x"^;e)\< Y. \V{x"^\^\x"^;e)\ 

Under the assumption [Mod] and from Corollary 3, we have 

E \V{x"'\^\x"';e)\\ =|Ao|£;(|y(O^^|$;0)|e-^Kl*-^*)) <+oo 

This means that for all G C Aq, there exists a random variable Y G L^{Pg*) 
such that 1-^2, G I < Y- Thus, under the ssumption [Mod] and from Theorem 4 
(Nguyen and Zessin (1979)) and from Corollary 3, we have Pe*— almost surely 

= £;(y(O^^|$;0)e-^Kl*^^*)). (6.3) 

We have the result by combining (6.2) and (6.3): Pe*— almost surely 

[/„(•; ^) ^ U{e) = £;(e-^Kl*^«) + V (0^^|$; 9) e-v(o"^l*;e*)) (5.4) 
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Step 2. £/„(<&;■) a contrast function 

Recall that Un{-',0) is a contrast function if there exists a function K{-,6*) 
(i.e. nonnegative function equal to zero if and only if = 9*) such that Pg* —almost 
surely Un{'f, 9) — Unif, 9) -^ K{9, 9*). From Step 1, wc have as n ^ +cx) 

K{9,9*) = £;re-^Kl*'«*)(e^Kl*-«)-^Kl*;«*) 
-(l + V (0^^|$; 9)-V (0^^|$; 9*) 

= £;(e-«*"''Kl^) (ef^-^*)"^""!*) - (1 + (0 - 9*fv (0^»))). 

(6.5) 

Let y G RP+^ \ {0}, and assume ifv {0"'\ip) = for A™ (g) Pe*-a.e. (m, ip). 
By assuming [Ident]. it follows that for i = 1, . . . ,£ (^ > p) i/'^v (0'"'|(^i) = 
for all {rrii, (pi) G Bi. From the injcctivity of the matrix with entries Wj(0™*|(/9i) 
for all {nil, ipi, . . . , mi, ipg) € Bi x ... x Bi it comes that y = which leads to a 
contradiction. Therefore, for 9^9*, the assertion {9 — 9*) v (0'"|(^) ^ holds 
for A™ Pe*— a.e. (m, </?). 

By noticing that the function t i-^ e' — (1 + 1) is nonnegative and is equal to 

zero if and only if i = 0, one concludes that the random variable e~ "'-" ' -' x 
/ (e-e*f -«(o^^|*) _^^^(.^ - 9*fv {Q^^\'^))\ is almost surely positive for 9^9* 

and equals to zero when 9 = 9*. 

Before ending this step, note that for any ip, Un{f] •) and K{-, 9*) are clearly 
continuous functions. 
Step 3. Modulus of continuity. 

The modulus of continuity of the contrast process defined for ail ip £Vt and 
all 7] > by 

W„(v?, ?7) = sup { |c/„(<^; 9) - [/„(^; 9') | : 0, 0' G 0, 1 10 - 0'| | < ??} 

is such that there exists a sequence {ek)k>i, with e^ ^ as A: — > +oo such that 
for all fc > 1 

P Aim sup [wnU,^>etS\=Q. (6.6) 

Let us start to write Wn (</?, \) < Wi^n {f, ^) + W2.n {f, ^) with 



W^l,ri ( "l^' T ) = sup ■ 



l^nl JA„ XM 

6',6»'g0, ||6»-0'|| < i 
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and 

W2,nU.^ = sup 1 1 Y. ^(a;">\^'";e)-^(^">\^'";e') 



■.e,e' e &,\\9-e'\\ <-}. 
kj 

From the assumption [Mod] it comes 

-.0,9' e 0,||0-6/'|| < J i 

Under the assumption [Mod], one can apply Theorem 4 (Nguyen and Zessin 
(1979)) to obtain as n — > +oo 

-^ / \\v{x^'^\^)Mdx"')^E{\\v{0''\^)\\), ioiPg.-a.e.^. 

So there exists nj, (fc) such that for all n > tiq (fc) we have 

for P^^-a.e. (^. By using the same arguments, one can prove that there exists 
njj ' (k) such that for all n > tIq [k) we have 

M^2,„ f^, 1) < I ^ (ll^ (o*» ||e-^Ki*^«*)) < ^ £; (||« (0^^|$) II) , 



k 

for Pg*— a.e. (^a. And so for all n > no(fc) = max(nQ (fc),nQ (fc)). we have 
Pe*— a.s. 

M/"™ f <^, ^ j < Wi.n iv^-j:) + ^2,n ( f^' -^ ) < ^' foi' ^e* - a-e. (f. 
with 5 = 4e^^ (llv (0*^1$) II). Since, 



hm sup <i W„ ((/?,- ) > - 



- nu{»-..(-l)^^} 

mGNn>m '^ V / > 

C U {w^n(^,^)>^},forP,.-a.e 



..e. ip. 
I \ «: / fi; I 

n>no{k) 

the expected result (6.6) is proved. 

Conclusion step. The Steps 1, 2 and 3 ensure the fact that we can apply a 

consistency result on minimum contrast estimators, see Guyon (1995). 
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6.3. Proof of Theorem 2 

Step 1. Asymptotic normality of Un {^;d*) 

The aim is to prove that for any fixed D, the following convergence in distri- 
bution as n ^ +00 

|A„|i/2 C/(i)(<i>;r) ^ AA(0,S(5,r)) (6.7) 



where the matrix S(D, 0*) is defined by (4.4). 

The idea is to apply to Un {^]0*) a central limit theorem obtained by 
Jensen and Klinsch (1994), Theorem 2.1. The following conditions have to be 
fulfilled to apply this result. For all j = 1, . . . ,p + 1 



(i) For aU i e i:\ E [ [LPV^l ($; 9*)] |$a- ) = 0. 



[a) For aUi c^Z^.E 



< +00. 



(Hi) The matrix Var (\An\^^^uii'\<^- 6*)) converges to the matrix S(D, 6*). 
Condition (i): From the stationarity of the process, it is sufficient to prove that 

E((LPL'-ll{<^;e*))\<^^.)=0. 



Recall that for any configuration ip 

(LPL^lliv-e*)) = -f ^,(x'"|(p)e-^(-'"l^^^*)M(d^") 

t;j(x™|(p\x™)v?(da;™). (6.8) 

AqxM 

Denote respectively by Gi{(p) and G2{(p) the first and the second right-hand 
term of (6.8) and hy Ei ^ E (G'i($)|(f>Ag = fA^)- From the definition of Gibbs 
point processes, 

E2^ ^ j , I Q(d^Ao)/^Ao(d^-'")lAo(:r)^^,(^™I^V")e-^('^^ol'^-o^"i. 

Denote by (p' = (lySAo, Va=)- Since Q is a Poisson process we can write 

E2 = „ ] , I Q{dv') j (^'(da;'-)lAo(x)z;,(x'"|(p\x™)e-^(^"ol'^As;e*) 
^Ao(</'A=) Jn Js 

^ ^ Q(V) 



ZAa{V>Al) Jn 

X / ^'(dx")lAo(x)i>,(x'"|^:^., U^AS \2;'")e-^(^^ol^"5^'*) 
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Now, from Campbell Theorem (applied to the Poisson measure Q) 
1 



E2 = 



^iidx"") 



Since from Slivnyak-Meeke Theorem, Q — Q;^, one can obtain 



AqxM 



1 



ZAoifA-JJn 

= -El 



AqxM 



-v(an'"lv;e*) -"^'(■^aoIvajj;^*) 



Condition (ii): For any bounded domain A and any configuration ip, one may 
write for j = 1, . . . , p + 1 



LPL^^i^-e*) 



< 4 



AxM 



Vj{x^^\ip)e' 



-^(-"l^'^*V(da;™) 



J2 t',(x">\x") 

From [Mod], both right-hand terms are integrable with respect to Pe*, which 
implies that for any domain A and in particular for A^, 



E 



(xPxW(cI>;r) 



< +CX). 



(1) 



Condition (in): let us start by noting that the vector LPL]^' {(p;0*) depends 
only on ip^^ for j £ IB (i, -S | . Let 



-(1) 



•(1) 



by using the stationarity of the process. From definitions, we can obtain 

Var(|A„|i/2C/W(<I>;0* 



|A„|-iVar(^iPi(,\'($;r) 



ije-fr. 



lAnl-^Ef E 



£;. 



■».j "^ E 



£;. 



«j 



,e/„nB(,|ej) ,e/„nB(.|ej)^ 
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E 



E LPL]l> ($; e*) LPL]l> ($; 6* 



(1) 

A, 



•(1) 



(1) 



= E[E[LPL)^>{^;9*)LPL'^>i<^;e*) |$a. 



'(1) 



-(1) 



= ElE [LPL'l> ($;r) |$A,j LPL'l' ($;r) 
= 
Denote by /„ the foUowing set 

We now obtain 

Var(|A„|i/2[7W(a>;r)) ^ \A^\-' J2 E ^'^ 



tei„ 



■ je/„nB iA & 



|A«rM E E 



E,,_ 



E^. 



*j 



+ E E 

Using the stationarity and the definition of the domain A„, one obtains 

|A"I"' E^ E ^^^- = \^n\-'\In\In\ ^ i^O. 

-^ ^{3,9*), 
as n — > +CXO and 



|A„ 






E,_ 



<|A. 



n M- 'H. 



E l^^ojI^O 



ieB 0, ^ 
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as n — > +00. Hence as n ^ +00 

Var(|A„|i/2C/W(cO;r)) ^ \K\-'Y1 E ^^^■ 

5-2 fcesf 0, & 



Step 2. Domination of Uj^'{^;6) in a neighborhood of 9* and convergence of 
Let j,k = 1, . . . ,p, recall that [U.nifj ^) ) is defined for any configuration 

Using tlie local stability property it conies 

^ ^ 3,k |A„| Ja„xM 

From [Mod], one can apply Theorem 4 (Nguyen and Zessin (1979)). It follows 
that there exists no such that for all n > uq 






< 2 X e^E (\vj{0^'\^)vk{0^^\^)\) for Pe* - a.e. ip. 



Note that from Theorem 4 (Nguyen and Zessin (1979)), for all 9 (and in par- 
ticular 9 = 9*), U_n {■] S) converges almost surely as n ^ -l-oo towards U}'^'{9) 
which is the (p, p) matrix with elements 

(c/(2)(0)) ^^E (z;,(0™|<f>)t;fc(0'"|$)e-^(o'"l^-«)) , for j, k=l,...,p. 

Let us underline that U} ' (9) is a symmetric definite positive matrix. Indeed, 
it is a positive matrix since for all y G RP"""^ 



y'u}'H0)y = Y.y^EU{o'^\<^H{o'^\<^)e-^i'"^''-''^) 



j-k 



£;((2/^t;(0*»)%-^(°"l*^^)) >0 



Vk 



And it is definite since, for all y G W+'^ \ {0} from [Ident], y^v {0"^\(p) = for 
A™ (g) Pg*—a.e. (m, ip) imphes y = 0- 

Conclusion Step Under the assumptions [Mod] and [Ident], and using Steps 
1 and 2, one can apply a classical result concerning asymptotic normality for 
minimum contrast estimators, see Guyon (1995) in order to obtain (4.3). 
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Now, the result (4.5) is proved in three substeps: 

(i) We first prove that the matrix S(0*) = S(Z3, 0*) is a symmetric definite 
positive matrix. From Equation (6.9), it is sufficient to prove that the matrix 
Yai{\Kn{D)\~^/'^LPL^ ' — ($; 6*)) is definite positive for n large enough. Let 
y ^MP \ {0}, the aim is to prove that 

V := j/^Var {\K{D)\-^'^LPL^ll-^ ($; 0*)) y > 0. 

Let A = U. / r^i\ Ai(D), using the same argument of Jensen and Kiinsch 
(1994) (Equation (3.2)), one can write 



V > \An{D)r' ^ ( Var ( y^LPL^l^-^ ($; r ) |<i>A,(l5),M (2 



1)Z^ 



Note that for i y^ j G In 



Cov y'LPL^']- (^; r ) , y'LPL^l\- (^; 6*) \<^.,T^,.k i (2 



(1) 



^,(D 



A,{D) 



-Afc(D)' 



1)Z^ 



= 



due to the independence of LPL\'—^ (^;0*) and LPL^,' —^ (^',d*) condi- 



tionally on ^A.fm''^ ^ (2 = + 1)^^ when ij e /„ n (2 £ + 1)2^ and 



-(1) 



LPL ' ,—^ (ip:6*) or LPL ' ,—^ (lo-.O*) is constant when either i or j 4 In H 



(2 



1)Z^. As a direct consequence, 



V > |A„(Z?)|-iiS Var /^LPi^i)_j($;0*)|$^^^-j,fc^(2 
= |A40)r^ ^i? (Var (/iPi^^_^_^ ($;r) |$^^(ij),M (2 






+ 1)Z^ 



+ 1)Z^ 



V V 'e^"nB(.,[g]) 

+ |A„(0)r^ J2 

£ I Var I / ^ LPL« _ ($;0*) |$,,(^,,fc ^ (2 



+ 1)Z^ 



+ 1)Z^ 
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Following the proof of Step 1, condition {Hi) one may prove that the second 
right-hand term tends to as n — > +cxd. Therefore by using the stationarity, we 
have for n large enough 



V > 



r|An(D)r 



/nn(2 






+ 1)2 



X E (Var {fhPL'^S' {^-9*) |$^^(o), 1 < |fc| < 2 



D-^ |/n n (2 



+ 1)Z^| 



X i; (Var (/LPLW ($;r) |<&^^(7j),l < |fc| < 2 



D" 



+ 1 



X i; (var (fLPL^^^ ($;r) |$^^(o), 1 < |fc| < 2 

Assume there exists some positive constant c such that Pg*—&.s. 

•(1) 



are (for 



ifLPL^' ($;0*) = c when the variables $Afc("D)'l < 1^1 < 2 
example) fixed to 0. By assuming [SDP] it follows that for any ipi G Ai for 
i ^ 0,...,e (with i >p),y^ (^^^T ^'^*' ^*) " ^^^T ^"^O' ^"^0 = °- ^'""^^ 
for all ((/3o, • • • , 'Pi) G ^0 X ■ • ■ X ^f- the matrix with entries ( LPLj- [ipi] 6*) ) — 

(LPL-— {ipo;6*)] is assumed to be injective, this leads to y = and hence 



to some contradiction. Therefore, when the variables <I>^ (-q), ^ < \k\ < 2 

are fixed to 0, the variable y'^ LPLi- ($; 0*) is almost surely not a constant. 
Hence, S(0*) is a symmetric definite positive matrix. 

(m) Convergence of ^^-^{ip; D^ , D , 6) . 

Let us recall that for any ip ^ il, D^ > D and G we define 



sj^ji^^A^) 



We also define 

X,{ip) := X.iip)''^' 



je/„nB 



E {LPL^^^{^;e))^[LPL^l]{^;e 

(4fl) 



yi(¥') := X,((^) - E{X,{<^)) and y„(v?) = |/„r' E.g/,. >'»(v)- Since one may 
notice that ^(Xj($)) = i52 (s(i5, 6)) , we have 

V / k,£ 

Y^{^) = & {%^{^-D\b,e)~UD,e) 
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Thus, the aim is to prove that as n ^ +00, Yn{ip) -^ for Pg* — a.e. tp. Since 
the process {Yi,i G Z^} is stationary, it is sufficient to prove, see e.g. Guyon 
(1995) 

• (a) E{Yo{^^) <+oo 

. (6)£;(|/„|F„($)2) <+w. 

(a) We leave the reader to verify that [Mod] ensures this integrabihty con- 
dition. 

. Hence, 

D 



(6) Note that Yi{ip) depends only on ip/\j for j G B f i 



D^ 



D 

by choosing j G /„ fl IB (j,aDv £))^ with ao^^.D — old"^ ,d{D) := 2 
the variables Yi and Yj are independent. Then, we obtain 

E{\In\Yn{<^f) = ^ E E{Y,{^)Y,{<^)) 

= toE( E E{YmY,m) 

+ J2 E{Y,{^)Y,{^))] 

= ^E E E{Ymy,m) 

E £;(ro($)n($))<(2az,v,B + i)£;(ro($)'). 

fee]B(o,Qov.B) 

Therefore, for aU Z)^ > I? and for all G 0, we have for Pg. — a.e. (^ as n ^ +00 

S„(^;I?^5,0)-^S(5,0) = S(0). (6.10) 

(iii) Since for any (^, the functions C/„ (<p; ■) and S„((y9; D^, l), •) are con- 
tinuous, it follows from Step 2 and (6.10) that one obtains for Pe*— a.e. (^, as 

n — > +CXJ 

C/(f)(<^.0)^[7(2)(0*) and S„((^;i?^,i5,0)^S(r). 

Finally, note that the previous convergence also implies that for n large enough 
S„(<1>; Z?^, Z?, &) is almost surely a symmetric definite positive matrix. 
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